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A formally exact expression for the tunneling current, for its separation into specular and diffuse 
components, and for its directionality, is given for a thick tunnel junction with rough interfaces in 
terms of the properties of appropriately defined scattering amplitudes. An approximate evaluation 
yields the relative magnitudes of the specular and diffuse components, and the angular dependence 
of the diffuse component, in terms of certain statistical properties of the junction interfaces. 

PACS numbers: 73.40 Gk, 74.50. +r 
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INTRODUCTION The study of the quantum me- 
chanical tunneling of electrons between two metallic elec- 
trodes separated by a thin barrier is an important method 
for investigating condensed matter systems (e.g. see Ref. 
0). Although the vast majority of tunneling experi- 
ments have been carried out on tunnel junctions whose 
interfaces have a significant roughness, the impressive 
theoretical literature |Q treating the properties of differ- 
ent types of tunnel barriers and tunneling mechanisms 
has almost without exception (see however Ref. 0j) dis- 
cussed only the case of flat tunnel junctions. This article 
presents the first detailed theory of tunneling appropriate 
for tunnel junctions with rough interfaces. The potential 
significance of such a development is apparent from one 
of our conclusions, namely that for junctions where the 
interface roughness fluctuations exceed an electron wave- 
length in magnitude, the contribution of the diffuse trans- 
mission of electrons to the tunneling current dominates 
the specular transmission that is usually calculated. 

A central idea in the flat interface theory of tunneling 
is that for thick barriers the electrons which dominate the 
tunneling are those whose momenta are directed close to 
the forward direction -|||. This "tunneling cone" ef- 
fect is the basis for attempts to determine the anisotropy 
of the superconducting energy gap (see page 126 of Ref. 
JjJ), and has also recently been invoked in the explanation 
of tunneling phenomena high temperature superconduc- 
tors [^]-|| where the spectrum of quasiparticle excitation 
energies is highly anisotropic. The investigation carried 
out below of tunneling directionality in the case of rough 
interfaces (where flat interface tunneling cone ideas are 
not applicable) thus has important implications for these 
studies. 

The theory of wave scattering at rough surfaces is a 
highly developed subject || with applications in many ar- 
eas of physics. Below, some established ideas from these 
studies, such as the use of certain scattering amplitudes 
and of ensemble averages over the random variables de- 
scribing the rough interfaces, are used to derive a formal 
expression for the tunneling current and to separate it 
into specular and diffuse components. This expression is 



then evaluated within the framework of two complemen- 
tary classical approximation schemes, a small perturba- 
tion method valid for roughness fluctuations smaller than 
the electron wavelength, and a quasiclassical approxima- 
tion (implemented via the tangent plane method) valid 
in the opposite limit. The approach of this article is thus 
quite different from a previous discussion of diffuse scat- 
tering in tunneling |Q| which has no way to separate the 
specular from the diffuse scattering, to calculate their rel- 
ative magnitudes, or to investigate the factors influencing 
directionality in the case of rough interfaces. 

The SUMMARY AND CONCLUSIONS section at the 
end of the paper gives an overview of the main results. 

FLAT TUNNEL JUNCTION INTERFACES AND 
THE TUNNELING CONE Consider an electron tun- 
neling from one metal to another through an insulating 
barrier. In the prototypical problem the electron is 
described by the Schrodinger equation 

(-7i 2 /2m)VV + V(z)tp = Eip. (1) 

The potential V(z) is shown in Fig. [ij The energy of the 
electron E lies between and Vb so that the insulating 
slab is a classically forbidden region. 
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FIG. 1. Barrier potential for an MIM' (Metal Insulator 
M'etal) tunnel junction. The random height functions as- 
sociated with interfaces at z — — t and z — are hi(r) and 
/i2(r), respectively. 

The electron wave function in the insulating region has 
the form 



V>/(R) = e ikr (Be- RlkZ + B'e KlkZ ), (2) 
K Jk = [nj + fc 2 ] 1 / 2 , kj = [(2m/h 2 )(V - E)] 1 ' 2 . (3) 



Throughout this article three-dimensional vectors are de- 
noted by boldface uppercase letters and two-dimensional 
vectors by boldface lowercase letters. Thus, R = (r, z). 
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The current tunneling across a thick (njt 3> 1) barrier 
in the presence of an applied voltage V is 

J z = || dE[f(E) - f(E - eV)] J -0^D(E, k), (4) 

D(E,k) = fl (k)ezp(-2/sjt)27r(Afci) 2 P Afc ,(k), (5) 
P Afcj (k) = ^(Akjfr'expi-ky^Ak!) 2 }}, (6) 

where D is the transmission coefficient and Pa/c 7 is a nor- 
malized (i.e. J PAfcj(k)cPfc = 1) two dimensional Gaus- 
sian function of width Aki = [kj / '{it)] 1 / 2 . The factor 
exp(—2Kjt) in Eq. |5| comes from the exponential decay 
of the wave function in the insulating layer, whereas the 
Gaussian function reflects the fact that for thick barriers, 
only electrons which have momenta close to the forward 
direction contribute to the tunneling current. (This last 
property is often called the "tunneling cone" effect.) The 
prefactor g is a relatively weakly varying function of mo- 
mentum and energy of order of magnitude unity, and can 
usually be neglected in calculations of the tunneling cur- 
rent (e.g. see the discussion on p. 21 of Ref. §). Below, 
g and its analogues will be put equal to unity. 

ROUGH TUNNEL JUNCTION INTERFACES - 
FORMAL THEORY Consider an electron which is in- 
cident on the tunnel junction from the metal M in Fig. 
[l]. The electron passes into the classically forbidden re- 
gion at z = — t and then, after having been attenuated 
by a factor exp(-nikt), it arrives at the z = interface 
where it is finally transmitted into the metal M' . Be- 
cause Kik is given by Eq. ^[ independently of whether 
the interfaces are rough or not, it is primarily electrons 
with small parallel momentum components k that arrive 
at the interface z = 0. If the interface at z = is rough, 
however, it will impart a random parallel component of 
momentum to electrons entering the metal M'. Thus 
the tunneling electrons will sample states in M' having a 
wide distribution of parallel momentum components, in 
spite of the tunneling cone effect in the insulator. The 
essential problem is to find how the transmitted current 
from electrons of wave vector k in the metal M is dis- 
tributed over the various wave vectors k' in metal M'. 

The calculations below use the methods described in 
Rcf. [|| , adapted here to the problem of tunneling. Con- 
sider first a general description of the transmission and 
reflection of a plane wave incident from medium 1 onto 
a rough interface separating media 1 and 2. The wave 
reflected back into medium 1 is written as a linear com- 
bination of waves with all possible parallel momentum 
components, as is the wave transmitted into the medium 
2. The wave functions in media 1 and 2 are thus 

WR) = ^T+ / T^u^k)-^, (7) 
Ilk J v ' Ilk' 

/ fly p'KjR 
^(k'.k)^-. (8) 



Here K tt = (k, q a k), a — 1,2 where q a k = inik if a refers 
to the insulating region, and q a k — (Kfj — k 2 ) 1 / 2 with 
Km = [(2m/h 2 )E] 1 / 2 if a refers to a metallic region. The 
division by q 1 ! 2 in Eqs. [7] and || represents the conven- 
tional normalization of the plane waves The quanti- 
ties Sii and S\2 are called scattering amplitudes. 

For the tunneling problem, the basic scattering am- 
plitudes are Si t M (i-e. 1 = M and 2 = / in Eq. ||) and 
Sm'J describing the transmission of an electron from the 
metal M to the insulator /, and from the insulator I to 
the metal M' , respectively. An important simplification 
in the calculation of S^m is that for the thick junctions 
considered here, only the exponentially decaying waves 
need be considered in the insulating region, and the ex- 
ponentially increasing waves can be neglected (e.g. see 
Rcf. |0|). It can be shown that in this approximation the 
scattering amplitude Sm'.m describing the transmission 
from metal M to metal M' is given by 

J d 2 fc"5M',z(k',k")5 7 , M (k" ) k)e i ^"" -«"*>*. (9) 

The metal-insulator interfaces are given in terms of 
the random functions h a (r), a — 1,2 by the equations 
z = — t + hi(r) and z — ^ia(r). The ensemble average 
of each h a (r) is taken to be zero, so that the average 
interfaces are flat. The boundary conditions satisfied by 
the wave function are that both the wave function and 
its normal derivative are continuous on both interfaces. 

The ensemble average of Eq. || must give an average 
wave function ip2 corresponding to flat interfaces; the av- 
erage scattering amplitude thus has the form 52i(k', k) = 
V2i(k)<$(k' — k). The scattering amplitude is now writ- 
ten as the sum of its average value and a fluctuating part: 

5 2 i(k',k) =F 2 i(k)(S(k'-k) + AS 2 i(k',k). (10) 

Furthermore the correlation function of the scattering 
amplitude fluctuations can be written in the form 

AS 2 i(k",k)A5 2 i(k',k) = <r 2 i(k',k)5(k" - k'). (11) 

Given that the ensemble average of the current density 
normal to junction in the metal M' can be calculated us- 
ing the formula J z = (h/m)Im('ip M ,d'ipM' /dz), the trans- 
mission coefficient D(E, k) appearing in Eq. |] can now 
be found using Eqs. §, [L0| and [□], with the result that 

D(E,k) - \V M >M(k)\ 2 + J a M , M {k',k)d 2 k', (12) 

where the integration is restricted to |k'| < Km- From 
Eq. it is clear that the fraction of the incoming cur- 
rent in M with parallel wave vector k transmitted into 
states in d 2 k' is cfm'.m (k', k)c? 2 /c' whereas the fraction 
transmitted without change in the parallel component of 
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momentum is \V m' .M(k)\ 2 . In terms of quantities char- 
acterizing the two junction interfaces, one finds 



|^M',M(k)| 2 = |V M ',j(k)Vj,M(k)| 5 



(13) 



and 



a WM {k',k) = a M '/(k',k)|y / , M (k)| 2 e- 2K - t 
+ |F M M(k')| 2 <T JM (k',k) e - 2re «' i 

+ yd 2 fc"a M , 7 (k',k") e - 2K ""V /M (k",k) (14) 

The diffuse contribution to the tunneling current, Eq. 
|l4| , contains contributions in which the transmission is 
diffuse at one interface and specular at the other, as well 
as a contribution (the last term) which is diffuse at both 
interfaces. 

THE SMALL PERTURBATION METHOD The 
small perturbation method works when the flat surface 
problem (i.e. the problem for h a (r) — 0) is a good first 
approximation. The corrections are calculated by ex- 
panding in powers of h a (r) . The quantities h a (r) appear 
in the calculations because expressions such as Eqs. [?] 
and |^ are evaluated at the interfaces z = — t + hi(r) and 
z = /12(f) when applying the boundary conditions. Thus 
h a (r) appears in expressions such as exp[— ni k h a {v)] and 
exp[iqMkh a (r)], and expansions in powers of h a (r) will 
be expansions in powers of the parameters [Ki k h a (r)] and 
[gMkh a (r)]. The quantities qMk and Kj k are of the order 
of magnitude of 2tt/X where the electron's wavelength A 
is expected to be comparable in magnitude to the lattice 
constant. The small perturbation approach will therefore 
be valid only when the root mean square fluctuations in 
h a (r) are smaller than a lattice constant, i.e. for atomi- 
cally flat interfaces. Since the results of the section on flat 
interfaces are already a good first approximation when 
the small perturbation method is applicable, no further 
results of this approximation will be given. 

THE TANGENT PLANE APPROXIMATION This 
section evaluates the transmission coefficient D(E, k) oc- 
curring in Eq. || for the tunneling current within the 
framework of the tangent plane approximation Q . This 
approximation works best when the spatial scale of the 
roughness is larger than the electron wavelength, and is 
thus complementary to the small perturbation approach 
outlined in the previous section. 

Consider first the general case of the transmission of 
a plane wave from medium 1 to medium 2 across the 
random interface z = h(r), which is described in terms 
of Eqs. and The method begins with a mathemati- 
cal formulation of Huygens' principle in which the wave 
function of the electron in the medium 2 is given in terms 
of its value and the value of its normal derivative on the 
interface z = h(r), namely, 



dn' 



G (R' -K)dS'. 



(15) 



Here S 1 is the surface z — h(r), R' is on this surface, and 
d/dn' is a normal derivative into medium 2. Also, the 
Green's function Go(R' — R) satisfies the equation (V 2 + 
iff )Gq(R' — R) = <5(R' — R) and can be represented as 



G»(R) = - jp 



exp[-ik ■ r + ig 2 fc|.z|] 



d 2 k. 



q2k 



(16) 



The next step is to find the electron wave function, 
V> 2 (R'), at points z = h(r) on the interface. This is 
done in the tangent plane approximation by consider- 
ing a given point on the interface, constructing a tan- 
gent plane there, and then considering the reflection and 
transmission of the incoming plane wave (which is taken 
to be the first term in Eq. [?]) at this tangent plane. This 
gives ^(R') in terms of the amplitude and phase of the 
incoming plane wave, and this result can be combined 



with Eqs. H, and 16 to yield 



5 2 i(k',k) = 

J exp[-i(k' - k) • r + z{q 2k - q lk )h(r)]^. (17) 

Here, a complicated function of the wave vectors of order 
unity, and analogous to the prefactor g in Eq. |^, has been 
omitted. 

The quantities Vm',m and <tm>,m necessary for an 
evaluation of the transmission coefficient D(E, k) (Eq. 
[l2] ) can now be evaluated by combining Eqs. |[ [10|, |ll| 
and [l?]. In carrying out the necessary ensemble averages 
the function h(r) is assumed to be Gaussian, and the 
theorem exp(-h) = exp(—h 2 /2) (valid for any Gaussian 
variable having h = 0) is used. The results are 



IV 



M'.M 



(k)| 2 =F/3 1 /3 2 P Afc7 (k) 



(18) 



and 



<7 M ',M(k',k) = F/hPAinQs! ~ k)P Afcj (k) 
+ F[3 2 P Akl (k')P Akl (k' -k) 



F J d 2 k"P Ak2 (k' - k")P Akl (k")P Akl (k" - k). (19) 



where 



F = 2 7 r(Afc / ) 2 e- 2K '*e 2K '(' l ?+' l i), 



(Ak a ) 2 



e -{^+K 2 M )hl ^ 

(« 2 + K 2 M )& 



(20) 

(21) 
(22) 



The P Aka 's are the normalized Gaussian functions de- 
fined by Eq. p| w ith widths Ak a given by Eq. ^2] where 
s a = (dh a /dx) 2 — (dh a /dy) 2 is the mean square slope 
of the roughness. 
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SUMMARY AND CONCL USIONS The approach to 
tunneling theory introduced above allows a calculation 
of the consequences of rough tunnel junction interfaces 
on the tunneling current and on its directionality. The 
general formula for the tunneling current is given by Eqs. 
f|, [l^, HI and |l4|. Eq. [l^ shows the separation of the 
current into specular and diffuse parts, and the scattering 
cross section om'M gives the directional dependence of 
the diffuse part. Eqs. [l3|, and [l4| reduce the problem to 
the determination of the transmission properties of the 
individual junction interfaces. These expressions give a 
formally exact theory of tunneling for thick rough tunnel 
junctions, and can be evaluated using any appropriate 
approximation scheme. 

The small perturbation method treats the roughness as 
a perturbation of a flat interface model, and shows that 
flat interface models represent a good first approximation 
when the amplitude of the roughness fluctuations is less 
than the electron wavelength (which normally requires 
atomically flat interfaces). 

The results obtained in the tangent plane approxima- 
tion show that for rough tunnel junction interfaces (i.e. 
surface height fluctuations significantly greater than the 
electron wavelength) the transmitted current is nearly 
totally diffuse. To see this recall that the functions P oc- 
curring m Eqs. H and HI are normalized Gaussians. Thus 
the relative weights of the different contributions in Eqs. 
|j~8| and [l9] to the tunneling current are determined by the 
factors (3 a . This means that for root mean square fluc- 
tuations in the height functions h(r) much greater than 
an electron wavelength, i.e. such that the factors (3 a are 
small, the purely diffuse contribution, namely that last 
term in Eq. [l9| dominates. 

Now examine the directionality in the rough junction 
case where the tunneling current is dominated by the last 
term in Eq. For sufficiently thick tunnel junctions, the 
factor PAfcj(k") = 6(k") and the integration over k" is 
easily carried out. The incoming and outgoing electrons 
contributing to the tunneling current thus have their par- 
allel moment components within Ak\ and A&2 (see Eq. 
p2|) of zero, respectively. For root mean square (rms) 
roughness slopes s a which are of order unity or not too 
much smaller, there is no significant directionality of the 
tunneling. On the other hand for rms roughness slopes 
much less than unity, the smaller the roughness slope, 
the closer to the forward direction is the momentum of 
electrons contributing to the tunneling current, both for 
incoming and outgoing electrons. 

It is of interest to examine the physical reasons for 
the dominance of the diffuse component of the tunnel- 
ing current in the case of rough junctions. As for flat 
junctions, the tunneling current is reduced by the fac- 
tor exp(—2nit) depending exponentially on the average 
thickness t of the junction, (see Eq. |2^). This effect is 
reduced by the interface height fluctuations, which give 
regions where the potential barrier has a smaller than 



average thickness. Hence the factor exp(2n 2 I (h\ + h\)) 
in Eq. |2(]. The reduction of the attenuation due to bar- 
rier thickness fluctuations is not as great for the specular 
component of the transmission, as indicated by the fac- 
tor exp(— Kj/Iq) in (3 a . The other factor contributing to 
[3 a , exp(—Kffha), gives the reduction in the specularly 
transmitted component of transmission due to destruc- 
tive interference of waves with the different phase lags 
due to having traveled different distances in the insula- 
tor. These two factors combine to make the specular 
component of the tunneling negligible relative to the dif- 
fuse component for sufficiently rough interfaces. Clearly, 
tunneling theory must account for the roughness of tun- 
nel junctions in order to correctly describe the dominant 
diffuse contribution to the tunneling current. 

This article has given a formally exact expression for 
the tunneling current valid for thick, rough tunnel junc- 
tions, has shown that for rough tunnel junctions the dif- 
fuse component of the tunneling current dominates the 
specular component, and has also shown that even when 
the tunneling current is entirely diffuse, a tunneling cone 
effect can exist if the root mean square roughness slope 
is sufficiently small. 
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